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INTRODUCTION 

Under  Contract  Nonr  3672(00)  with  the  Office  of  Naval  Research, 
Cornell  Aeronautical  Laboratory  is  conducting  a  program  to  investigate 
precipitation  processes.  The  program  is  divided  into  two  phases  and 
includes:  (1)  studies  of  cloud  and  precipitation  processes  using  pulsed 
Doppler  radar,  and  (2)  measurements  of  atmospheric  space-charge  in 
the  vicinity  of  clouds.  This  is  a  summary  of  the  technical  activities  under 
item  1  above. during  the  first  half  of  the  program,  which  ended  31  January 
1963. 
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DOPPLER  RADAR  INVESTIGATIONS 


The  purpose  of  this  phase  of  the  project  is  to  use  a  Doppler  radar 
svstem  for  quantitative  measurements  of  characteristics  of  the  wind  field 
ar  d  the  hydrometeors  associated  with  precipitation  processes  in  the 
atmosphere.  It  was  demonstrated  in  an  earlier  project  (Rogers  and  Piliff, 

1962)  that  radar  techniques  could  be  used  for  measuring  drop-size 
distribution  in  rainfall  under  stable  atmospheric  conditions.  It  is  intended 
to  extend  these  techniques  in  the  present  project  to  provide  measurements 
of  the  growth  rate  of  rain  drops  that  occur  in  relatively  stable  atmospheric 
conditions  and  to  make  simultaneous  measurements  of  drop-size  distributions 
and  the  speed  of  updrafts  in  convective  clouds*.  Further  we  plan  to  investi¬ 
gate  the  relationship  between  turbulence  and  the  precipitation  process. 

At  the  time  this  program  was  initiated  the  relatively  insensitive 
Doppler  radar  used  for  the  Mt.  Withington  experiments  (Rogers  and  Pilie, 

1962)  was  planned  for  use  on  the  project.  Shortly  thereafter  the  Laboratory 
provided  funds  for  extensive  modification  and  improvement  of  the  radar. 

The  modification  program  encountered  difficulties  which  as  indicated  by 
Section  B  oC  this  report  are  only  now  being  ironed  out.  At  the  present  time 
we  are  confident  that  the  radar  will  be  in  operation  in  the  field  by  approxi¬ 
mately  15  March. 

During  the  time  the  equipment  has  been  inoperative  we  have  completed 
the  analysis  of  certain  problems  which  arise  in  interpreting  Doppler  radar 
data  from  weather  targets.  This  analytical  work  has  led  to  two  papers  which 
are  attached  as  appendices.  One  of  the  papers  will  be  published  in  the  March 
issue  of  the  Journal  of  the  Atmospheric  Sciences  and  the  other  was  presented 
last  October  at  the  University  of  Michigan's  Second  Conference  on  Remote 
Sensing  of  Environment.  The  conclusions  of  both  these  papers  are  summarized 


Since  the  inception  of  this  program,  a  report  has  been  given  of  such 
simultaneous  measurements  of  updraft  intensity  and  size  distributions. 
(Battan,  L.  J.  ,  et  al.  •  "Observations  of  showers  by  means  of  pulse 
Doppler  radar,  "  paper  delivered  at  211th  Meeting  of  Amer.  Meteor.  Soc., 
Jan.  1963,  New  York). 
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in  the  body  of  this  report  in  addition  to  an  analysis  of  the  interaction  between 
hydrometeors  and  the  wind,  which  is  essential  for  determining  how  accurately 
the  wind  can  be  measured  with  Doppler  radar.  Finally,  a  brief  description 
will  be  given  of  some  of  the  characteristics  of  the  improved  Doppler  radar 
system-. 

A.  Analytical  Work 

1.  Meteorological  information  contained  in  radar  signal  fluctuations^ 

The  most  complete  information  obtainable  from  fluctuations  in  a 
weather  radar  signal  is  the  distribution,  within  the  sampled  volume,  of  the 
radial  component  of  the  velocities  of  tht  hydrometeors  responsible  for 
scattering  the  electromagnetic  energy.  Coherent  Doppler  radar  is  needed 
to  obtain  this  velocity  distribution,  from  which  properties  of  the  target  such 
as  drop  sizes  and  wind  speed  can  be  inferred.  It  has  been  known  for  some 
time  that  less  detailed  information  can  be  derived  from  measurements  with 
ordinary  incoherent  pulsed  radars.  A  theory  for  interpreting  incoherent 
signals  was  worked  out  by  Fleisher  (1953)  on  the  basis  of  an  analysis  by 
A.  J.  F.  Siegert,  which  is  reported  by  Lawson  and  Uhlenbeck  (1950). 

Siegert's  analysis  can  also  be  used  to  explain  the  relationship  between  the 
power  spectrum  of  the  signal  received  by  a  Doppler  radar  and  the  velocities 
of  the  scatterers.  This  relationship  then  provides  a  basis  for  explaining 
coherent  or  incoherent  measurements. 

The  fundamental  relationship  between  the  measured  power  spectrum 
of  the  received  signal  and  the  radial  velocity  distribution  S (v.)  of 
the  scatterers  is  found  to  be 

Fjf)‘  sU-  Q-) 

where  -f '  is  the  frequency  to  which  the  received  signal  may  be  heterodyned 
for  convenience  in  spectrum  analysis,  ^m(f)  and  therefore  the  velocity 
distribution  are  available  in  Doppler  radar  experiments. 


A  detailed  treatment  of  this  topic  is  given  in  Appendix  A. 
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If  incoherent  radar  is  used  instead  of  Doppler  radar,  the  measurable 
property  of  the  signal  is  its  power.  The  spectrum  of  the  power  fluctuations 
is  determined  by  the  velocity  distribution,  but  in  general  the  relationship  is 
not  unique  so  that  this  spectrum  can  not  be  used  to  obtain  the  velocity  distribu¬ 
tion.  Only  in  the  special  case  of  a  symmetrical  distribution  of  velocity  can 
measurements  of  signal  power  yield  the  distribution.  The  rate  at  which  the 
fluctuating  power  crosses  its  average  level,  however,  affords  a  simple 
measurement  of  the  spread  of  the  velocity  distribution,  regardless  of  its 
exact  shape.  If  the  rate  of  average  power  level  crossings  is  L  ,  then 

L.  -  S.  <£/a  <2> 

where  0^  is  the  standard  deviation  of  the  radial  component  of  velocity 
and  A  is  the  radar  wavelength. 


2.  Interpretation  of  the  velocity  distribution 


The  basic  piece  of  information  available  from  Doppler  radar 
measurements  is  the  velocity  distribution  S(u.)  of  the  hydrometeors  in  the 
sampled  volume.  This  velocity  distribution  is  determined  by  the  wind  in 
the  sampled  volume,  the  response  of  the  hydrometeors  to  the  wind,  and 
their  motion  relative  to  the  air  due  to  gravitational  settling.  In  the  case  of 
vertical  viewing  in  rain  ar.d  in  the  absence  of  significant  updrafts,  the 
velocity  distribution  is  caused  entirely  by  gfavitational  settling.  Then  the 
drop-size  distribution  n(D)  is  related  to  the  measured  velocity  distribution 
by 

„/n\  S(u)  dU  . 

m  x  -tfr  Hv  ’  (3) 

where  u(D)  is  the  terminal  fall  velocity  of  a  drop  having  diameter  D 
(cf.  Rogers  and  Pilif,  1962). 
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When  turbulence  and  updrafts  are  present,  the  complete  drop- size 

distribution  cannot  be  computed  from  the  spectrum.  In  this  situation  it  is 

possible,  however,  to  obtain  an  estimate  of  both  the  speed  of  the  updrafts 

1 

and  the  size  of  the  largest  drops  present  .  The  low-velocity  end  of  the 
spectrum  is  attributed  to  the  very  small  drops  which  move  essentially  with 
the  wind.  The  high-velocity  end  is  attributed  to  the  large  drops.  Since  the 
vertical  speed  of  the  wind  and  the  downward  speed  of  the  large  drops  toward 
the  radar  are  both  measured,  the  size  of  the  drops  can  be  computed  assuming 
they  fall  with  terminal  velocity  relative  to  the  air. 

When  the  radar  is  pointed  horizontally,  the  velocity  component  due  to 
gravitation  is  not  observed,  and  all  the  velocities  in  the  distribution  can  be 
explained  in  terms  of  the  horizontal  wind  and  the  response  of  the  hydrometeors 
to  this  wind.  Likewise,  for  any  viewing  angle,  such  weather  targets  as 
clouds  and  snow  have  velocities  containing  only  small  contributions  from 
gravitation.  In  order  to  derive  information  about  the  wind  from  the  Doppler 
spectra  from  these  targets  requires  understanding  the  way  the  hydrometeors 
respond  to  the  wind. 

3.  .Response  of  hydrometeors  to  wind 

If  raindrops  and  snowflakes  move  exactly  with  the  wind,  then  their 

Doppler  spectra  represent  the  distribution  of  the  radial  component  of  wind 

in  the  sampled  volume.  The  mean  velocity  in  the  spectrum  would  correspond 

to  the  mean  wind  in  the  volume  and  the  spread  of  the  spectrum  would  provide 

a  measure  of  turbulence  and  shear  within  the  sampled  volume.  This  interpreta- 

3 

tion  of  the  mean  velocity  was  reported  by  us  previously  and  has  also  oeen 
used  by  Lhermitte  and  Atlas  (1961).  A  small  error  could  be  introduced  into 
these  mean  wind  measurements  if  strong  shear  i'l  present.  In  that  case, 
a  raindrop  fhlling  into  the  sampled  region  would  not  immediately  assume  the 
velocity  of  the  air  in  its  environment,  but  would  fall  some  distance  before 
establishing  equilibrium.  These  effects  in  measurements  of  mean  velocity 
would  usually  be  quite  small. 


5 


roject  GUSH  Occasional  Report,  Nov.  1961. 
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A  more  difficult  problem  is  to  describe  the  manner  in  which 
quantitative  descriptions  of  turbulence  can  be  obtained  from  Doppler  measure¬ 
ments.  For  this  application  it  is  necessary  to  know  how  well  hydrometeors 
respond  to  the  rapid  turbulent  fluctuations  in  wind.  The  forces  acting  on  a 
bydrometeor  are  assumed  to  be  gravitation,  buoyancy,  and  drag.  Of  these, 
only  the  non-linear  drag  force  is  troublesome.  The  drag  is  given  by 
J-CdAPc L  V*  where  Cp  is  the  drag  coefficient,  <«,  is  the  density  of 
ambient  air,  A  is  the  cross  sectional  area  of  the  hydrometeor  and  Vf 
is  the  relative  velocity  between  hydrometeor  and  air.  Two  kinds  of  weather 
targets,  which  in  a  sense  represent  limiting  cases,  can  easily  be  treated 
because  their  pertinent  equations  of  motion  can  be  linearized. 

The  first  of  these  special  cases  is  a  target  consiting  of  cloud  droplets. 
When  the  Reynolds  number  is  less  than  approximately  10,  the  drag  coefficient 
for  a  sphere  varies  inversely  with  velocity.  In  this  region  of  flow,  the  drag 
force  on  a  sphere  of  diameter  D  is  equal  to  3lT//-VrD  where  Jl  is  the  viscosity 
of  air.  Because  of  the  limitation  on  Reynolds  number  this  form  for  the  drag 
applies  only  to  droplets  smaller  than  about  0,1  mm  in  diameter,  which  includes 
almost  all  cloud  droplets.  The  equation  of  motion  for  one  of  these  small  drops 
is 

JTf  -*  -*  —* 

31T/JiVVr  -  ntf  i-  <f- 

In  order,  the  terms  on  the  right  are  the  drag  force,  gravitational  force,  and 

buoyancy  force.  Since  the  mass  of  air  displaced  by  the  droplet  is 

-3 

approximately  10  times  the  mass  W  of  the  droplet,  the  buoyancy  term  can 
be  neglected.  Then  the  equation  reduces  to 


where  -A, 9  SW^JL  J) / fll  .  Substituting  Vr*  ,  where  is  the  velocity 

of  ambient  air, yields 

+  Jk,  7  =  JL,  7„  -  J. 


£7 

dt 
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The  equation  for  vertical  components  is 

wit'  obvious  notation.  In  the  case  of  stagnant  air  0  •  then  the  steady 

state  solution  is  simply 


>4  *-£ 


the  terminal  fall  velocity. 


The  general  steady  state  solution  for  the  case  in  which  the  air  oscillates 
sinusoidally,  according  to  -ur'  €,  ,  is 


JLrZH 

Jk'+Tul 


Lost 

e 


This  result  can  be  expressed  in  terms  of  filter  theory:  the  input  or  driving 
function  is  W  while  the  output  is  the  drop's  response  V^.  .  The  ratio 
of  output  (adjusted  by  the  terminal  fall  velocity)  to  input  is 


W 


!-*■«- -j - 


^  (to), 


which  is  identical  to  the  response  function  for  a  low  pass  filter.  The  analogy 
can  be  carried  further  with  useful  results  by  computing  the  system  function, 
which  by  definition  is 

H(<o)  5  J2[u>) 


For  this  case 


(4) 


The  usefulness  of  the  system  function  lies  in  its  property  that  the  spectrum 
of  the  output  is  the  spectrum  of  the  input  multiplied  by 


The  system  function  corresponding  to  horizontal  motions  can  be 
shown  to  be  identical  to  (4).  Substituting  for  ,  and  setting  ^t=1.7  x  10 
in  cgs  units  leads  to 


Hfw)  - 


I _ 


D  iVi  cm , 
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as  the  system  function  in  terms  of  droplet  size,  for  diameters  smaller  than 
0. 1  mm.  Figure  1  showu  this  function  for  a  droplet  of  0. 1  mm  diameter 
which  corresponds  to  the  largest  cloud  droplets  and  the  largest  permissable 
size  for  this  approximation.  It  is  seen  in  the  figure  that  frequencies  as  high 
as  one  cycle  per  second  are  followed  to  within  95%  by  the  droplet. 


A  second  special  case  for  which  the  equation  of  motion  can  be 
linearized  is  applicable  to  the  vertical  component  of  a  drop  having  high 
terminal  fall  velocity, and  was  first  recognized  by  Stackpole  (1961). 


The  general  equation  of  vertical  motion  for  raindrops  having  Reynolds 
numbers  higher  than  about  300  is 

Av. 


=  AJVrtVr 


where  -  CpPoJf V/gt*  and  KL*—  l4-  •  This  form  is  appropriate  for  drops 
larger  than  1  mm.  It  is  assumed  that  the  main  contribution  to  Vf  is  the 
terminal  fall  velocity  ^  .  Turbulent  deviations  are  then  taken  to  be  perturba¬ 
tion  quantities,  negligible  in  higher  orders.  The  system  function  is  then 
found  to  bo  . 


It  is  more  convenient  to  have  the  system  function  in  terms  of  drop  diameter 
rather  than  terminal  velocity.  The  data  of  Gunn  and  Kinzer  (1949)  can  be 
used  for  this  conversion.  The  system  function  for  various  drops  is  shown 
in  Figure  2. 


One  concludes  from  the  results  shown  in  Figs.  1  and  2  that  all  cloud 
droplets  respond  satisfactorily  to  be  considered  good  tracers  for  gusts,  but 
that  large  raindrops  do  not.  For  example,  the  5  mm  drop  is  seen  to  have  a 
spectral  response  of  only  10%  to  a  1  cps  driving  force,  a  value  entirely 
inadequate  for  a  tracer  material. 
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Still  unaccounted  for  are  the  drops  too  large  for  the  Stokes  flow 
approximation  and  too  small  for  the  perturbation  approximation.  These  are 
drops  in  the  approximate  interval  0.1-1  mm.  The  equation  of  motion  must 
be  solved  in  its  non-linear  form  in  this  region.  Further,  since  the  equation 
in  non-linear,  the  concept  of  a  system  function  does  not  apply.  Numerical 
solutions  to  the  equation  carried  out  on  an  analog  computer  indicate  that  only 
the  smallest  of  the  drops  in  the  0.1-1  mm  range  can  be  considered  adequate 
tracers  for  moderate  wind  gusts. 

This  same  non-linear  analysis  must  be  applied  to  snowflakes,  as  their 
equations  of  motion  do  not  lend  themselves  to  linearization.  The  conclusion 
of  the  analysis  is  that  dry  snowflakes  are  satisfactory  tracers  for  gusts,  but 
that  large  wet  snowflakes  respond  too  slowly  to  give  a  good  representation 
of  the  wind. 

In  conclusion,  the  analysis  of  hydrometeor  response  has  shown  that 
targets  consisting  of  dry  snowflakes  or  cloud  particles  are  nearly  perfect 
indicators  of  all  scales  of  atmospheric  turbulence.  Wet  snow  and  rain, 
however,  respond  poorly  to  high  frequency  gusts  and  can  only  be  used  as 
indicators  of  scales  of  motion  characterized  by  low  frequencies,  as  for 
example, the  "mean  wind"  of  scales  comparable  to  those  sensed  by  a  radio¬ 
sonde  balloon. 

4 

4.  Measurement  of  turbulent  energy 

For  those  weather  targets  consisting  of  particles  that  respond 
faithfully  to  the  wind,  the  velocity  distribution  of  the  particles  can  be  inter¬ 
preted  as  that  of  the  wind.  Therefore,  Doppler  radar  can  be  considered 
as  a  volume- sampling  anemometer.  \s  well  as  providing  an  estimate  of  the 
mean  wind,  the  Doppler  measurements  can  be  analyzed  to  obtain  the  total 
energy  in  the  radial  component  of  the  turbulent  wind.  It  is  required  that  two 
properties  of  the  Doppler  spectrum  be  recorded  continually:  the  mean  frequency 
(or  mean  radial  velocity)  and  the  variance  of  the  frequency,  as  indicated  by 
the  Doppler  spe'-trum.  These  quantities  are  denoted  re  spectively  by  and  (T*  . 


This  subject  is  treated  in  more  detail  in  Appendix  B. 
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Both  of  these  quantities  can  be  measured  in  real  time  using  analog  techniques. 
The  average  velocity  is  an  indicator  of  relatively  large  scales  of  turbulence, 
while  the  variance  is  an  indicator  of  the  relatively  small  scales.  The  total 
energy  in  all  scales  of  turbulence  can  be  obtained  by  analyzing  both  of  these 
quantities.  In  particular,  let  denote  the  variance  of  ( V )  over  several 

minutes  of  record;  and  let  denote  the  average  value  of  <j\ 

over  the  same  period.  Then  the  energy  E,  per  unit  mass  of  air,  in  the 
measured  component  V  of  turbulence  is  related  to  these  quantities  by 

+  <]•  (5) 

B.  Radar  Equipment 
1.  Introduction 

A  detailed  investigation  was  made  of  the  CAL  Doppler  radar  at  the 
beginning  of  this  project"*.  This  investigation  indicated  that  several  modifi¬ 
cations  and  additions  to  the  system  would  be  desirable  for  meteorological 
observations  such  as  planned  on  this  program.  Briefly  these  include: 

a.  The  installation  of  a  high  quality  stalo  in  order  to  improve  the 
short-term  frequency  stability  of  the  equipment  and  hence  the  Doppler 
velocity  resolution  capability, 

b.  The  installation  of  stablized  amplifiers  and  bias  power  supplies 
in  order  to  provide  accurate  and  stable  echo  cross-section  measuring 
capability, 

c.  Improved  transmitter  pulse  shape  and  peak  power  output, 

d.  Elimination  of  the  reference  signal  in  the  Doppler  spectra  by 
redesign  of  the  receiver, 

e.  Improvement  of  receiver  dynamic  range  and  linearity  characteristics. 


This  investigation  is  documented  in  a  CAL  memorandum:  B.R.  Tripp, 
MEMORANDUM  FOR  THE  RECORD,  "Summary  Report  of  DOP  Radar  Investi¬ 
gation  and  Proposed  Improvements,  "  13  August  1962,  BRT:pwm-26. 
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f.  Addition  of  an  average  frequency  tracker, 

g.  Elimination  of  spectra  aliasing  problems  by  altering  the  bandpass 
characteristics  of  the  receiver,  and 

h.  Other  more  minor  but  important  modifications  to  correct  drift  and 
spurious  responses  in  the  receiver. 

The  following  sections  discuss  briefly  major  aspects  of  the  radar 
system  and  related  instrumentation  now  available  for  this  project. 


2.  System  Design 

The  original  Doppler  radar  system  produced  a  strong  reference  signal 
component  in  all  spectra.  This  signal  was  caused  by  signal  leakage  from  the 
reference  oscillator  heterodyning  with  the  local  oscillator  signal.  Normally 
this  signal  component  is  of  no  concern.  However,  in  order  to  accurately 
measure  the  mean  frequency  of  Doppler  shifted  spectra,  it  was  necessary  to 
eliminate  or  suppress  this  leakage  components  to  below  the  receiver  noise 
level.  A  redesign  of  the  system  has  produced  an  acceptable  feedthrough  level 
on  the  basis  of  early  bench  tests  of  the  modified  receiver.  Final  performance 
will  be  measured  after  the  equipment  is  reinstalled  in  the  radar  van.  A  block 
diagram  of  the  modified  receiver  and  transmitter  (block  diagram)  is  given 
in  Figure  3. 

The  second  major  problem  was  the  short-term  frequency  stability 
of  the  receiver  stalo.  The  original  radar  used  a  free-running  reflex  klystron 
oscillator  that  was  unstable  and  hence  unsuitable  for  accurate  measurement 
of  Doppler  spectra.  This  problem  has  been  corrected  by  installing  a  stalo 

g 

with  a  short  term  frequency  stability  of  1  part  in  10  ,  or  100  cps  maximum 
drift  at  X-band.  During  time  intervals  corresponding  to  the  anticipated  signal 
propagation  time,  the  frequency  drift  should  be  one  order  or  more  of  magnitude 
less  than  the  maximum  permissible  value  of  100  cps.  An  optimistic  estimate 
of  Doppler  frequency  resolution  of  the  system  is  10  cps  for  target  ranges  of 
a  few  kilometers. 
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3.  Receiver  Characteristics 


Basically,  the  receiver  design  has  been  altered  to  reduce  spurious 
responses  and  to  extend  its  dynamic  range.  Several  very  narrow  bandwidth 
filters  have  been  added  prior  to  various  heterodyning  processes  to  minimize 
aliasing  errors  in  the  Doppler  spectra.  All  spectra  foldover  products  are 
down  by  at  least  55  db  in  the  Doppler  output  signal. 

A  typical  frequency  response  curve  of  the  modified  receiver  is  given 
in  Figure  4.  This  curve  shows  the  receiver  response  over  a  10  kc/s  total 
bandwidth.  The  slight  asymmetry  is  attributed  to  the  response  of  a  mechanical 
filter  which  is  included  in  the  system  in  order  to  achieve  the  rapid  fall-off 
of  frequency  response  outside  the  receiver  bandpass.  The  effect  of  asymmetry 
can  be  removed  if  necessary  from  the  Doppler  signal  by  multiplying  the 
Doppler  spectra  by  the  inverse  of  the  system  frequency  response. 

A  sample  calibration  curve  of  receiver  linearity  is  given  in  Figure  5. 
This  curve  illustrates  the  wide  dynamic  range  capability  of  the  modified 
Doppler  receiver. 

4.  Auxiliary  Instrumentation 

The  mean  frequency  of  the  Doppler  spectra  will  be  automatically 
measured  and  recorded  with  a  mean  frequency  tracker.  This  instrument  is 
currently  being  constructed  and  should  be  available  for  tests  within  the  next 
few  weeks. 

Other  instrumentation  which  have  been  added  to  the  radar  for  this 
program  include: 

a.  A  power  detector  and  pulse  integrator  for  measuring  target  cross 
section, 

b.  Detector  for  R-meter  tests, 

c.  Calibrated  readouts  on  all  system  controls,  and 

d.  Recording  facilities  for  making  records  of  all  data. 
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Figure  4  DOPPLER  RECEIVER  RESPONSE  -  12  Kc/s 
FILTER  POSITION 
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SIGNAL  INPUT  -  db» 


Figure  5  RECEIVER  LINEARITY  CURVE 
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ABSTRACT 


Statistical  properties  of  the  field  backscattered  from  a  weather 
radar  target  are  related  to  the  velocity  distribution  of  the  scatterers 
comprising  the  target.  The  complete  distribution  of  radial  velocity 
components  can  be  derived  from  the  received  radar  signal  using 
coherent  Doppler  techniques.  The  standard  deviation,  and  in  special 
cases  the  complete  form  of  this  distribution,  can  be  obtained  from 
incoherent  radar  measurements. 


The  material  presented  in  this  paper  was  developed  in  the  course  of 
several  programs  sponsored  by  the  Office  of  Naval  Research. 


A- 1 


VC-1606-P- 


1.  Introdu  .ion 


When  the  raindrops  or  snoflakes  comprising  a  weather  radar  target 
move  toward  or  away  from  the  transmitter,  the  backscattered  field,  observed 
at  a  point,  is  seen  to  fluctuate  randomly  in  both  amplitude  and  phase. 

It  is  well  known  that  these  fluctuations  can  be  analyzed  to  obtain  certain 
information  about  the  velocities  of  the  scatterers.  The  most  complete 
information  obtainable  is  the  distribution  of  the  radial  components  of 
the  velocities  of  the  scatterers.  To  measure  this  veloc:.,  distribution 
requires  the  use  of  carefully  designed  Doppler  radar  equipment.  Such 
measurements  have  been  irterpreted  in  terms  of  the  mean  wind  profile 
(Lhermitte  and  Atlas,  1961;  Lhermitte,  1962),  updrafts  in  showers 
(Probert-Jones  and  Harper,  1961),  and  raindrop-size  distributions  (Boyenval, 
I960;  Probert-Jones,  I960;  Rogers  and  Pilie,  1962). 

Before  Doppler  radars  became  available  for  meteorological  research, 
incoherent  radars  were  used  to  measure  some  properties  of  the  velocity 
distributions  in  weather  targets.  The  complete  velocity  distribution  cannot 
be  derived  from  these  measurements,  in  general,  because  the  phase 
information  is  lost  without  coherent  (Doppler)  techniques.  In  past  years, 
however,  considerable  attention  was  given  to  measuring  wind  shear  and 
turbulence  with  incoherent  radars;  these  efforts  are  reviewed  by  Battan  (1959). 
The  theoretical  foundation  for  interpreting  incoherent  ra<(ar  measurements 
was  worked  out  by  Fleisher  (1953)  on  the  basis  of  the  theory  of  A.  J.  F.  Siegert, 
which  is  reported  for  example  by  Lawson  and  Uhlenbeck  (1950). 

The  analysis  of  radar  signals  to  determine  the  information  about 
velocities  contained  therein  can  be  based  on  Siegert's  theory,  whether 
incoherent  or  coherent  radar  is  used.  In  view  of  the  increasing  use  of 
Doppler  radar  in  meteorological  work,  it  may  be  helpiul  to  show  how  the 
interpretation  of  Doppler  measurements  can  be  justified  and  to  show  the 
connection  between  Doppler  and  incoherent  measurements. 
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2.  The  autocorrelation  of  the  complex  amplitude 


The  backscattered  field  from  the  moving  collection  of  hydrometeors 
constituting  a  weather  radar  target  will,  in  general,  be  both  amplitude  - 
and  phase-modulated.  A  Cartesian  component  of  its  electric  vector  at 
a  fixed  point  in  space  can  therefore  be  represented  as 

E(t)  *  /\tt)  effete)  -  ,  (i) 

where  lOt-znf0  is  the  angular  frequency  of  the  carrier  wave,  and  Aft) 

and  <&(+)  are  real  functions  demonstrating  the  amplitude  and  phase  modulations. 

An  equivalent  representation  of  Elf)  is 

where  denotes  the  real  part.  This  representation  is  useful  because 
exponential  functions  are  often  easier  to  manipulate  in  calculations  than 
trigonometric  functions. 

A  third  representation,  and  the  one  we  will  use  to  describe  the  signal, 
is 

Eli)'  (2) 

The  quantity  6UJ  6  is  called  the  analytic  signal  associated  with  the 

real  signal  Eft)  .  The  complex  amplitude  €?tt)  is  given  by 

e(t)~  (t)  +  i  G  (+)  *  A(t)€  ) 

and  contains  both  the  amplitude  and  phase  modulation  effects.  It  is  made 
up  of  contributions  from  all  scatterers  comprising  the  target: 

4k  ’ 
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where  and  are  the  amplitude  and  phase  of  the  component  from  the 
M  scatterer.  This  complex  amplitude  is  not  a  measured  quantity,  but 
merely  a  mathematical  convenience.  Certain  properties  of  are  closely 

related  to  corresponding  properties  of  the  real  signal,  however,  and  computa¬ 
tions  are  more  readily  carried  out  on  the  analytic  signal. 

As  the  scatterers  move  about,  the  amplitudes  of  their  individual 
contributions  to  the  composite  field  are  assumed  to  be  constant  while  the 
phases  vary  according  to 


where  A  is  the  radar  wavelength  and  is  the  velocity  component  of  the  Jtfd 
scatterer  in  the  direction  of  the  transmitter.  Using  this  model  for  the 
scattered  field,  Siegert  derived  the  autocorrelation  function  &(t)  of  the 
complex  amplitude: 


&&)- 


eH)  _ 


(3) 


WJT c 

~0  -00 

where  overbars  denote  ensemble  averages  and  asterisks  denote  complex 
conjugates.  The  average  quantity  f*'  is  proportional  to  the  power  in  B(-f) 
and  is  defined  by 


U. 


V\/(f,U^  is  the  joint  probability  that  a  scatterer  picked  at  random  has  a 
velocity  component  toward  the  transmitter  in  the  interval  dv,  while  returning 
a  field  having  amplitude  in  df  .  G(t)  is  a  complex  function  of  'T  with  the 
demonstrable  property 


60 r)-  G*(-T)  * 
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this  property  is  a  sufficient  condition  for  the  power  spectrum  of  the  complex 
amplitude  to  be  real. 

Of  the  various  assumptions  used  by  Siegert  in  deriving  equation  (3), 
several  are  of  questionable  validity  when  applied  to  weather  targets.  Although 
it  is  not  the  intention  of  this  paper  to  assess  the  consequences  of  these 
assumptions,  it  should  be  mentioned  that  both  theoretical  and  empirical 
evidence  suggest  that  no  significant  errors  arise  from  applying  equation  (3) 
to  weather  signals. 

3.  The  power  spectrum  of  the  complex  amplitude 

The  available  information  about  the  velocities  of  the  scatterers  is 
contained  in  the  probability  Wty iu)  which  appears  in  equation  (3)  for  the 
autocorrelation  function.  This  probability  is  more  simply  related  to  the 
power  spectrum  of  the  complex  amplitude  than  to  its  autocorrelation.  The 
spectrum  f(-f)  is  related  to  the  autocorrelation  by  the  transform  pair 

F(f)-Jsfr)e'17r'fTc/T 

£(t)  =J  jit  . 

'oo 


The  two-sided  Fourier  transform  is  used  because  (rtf)  ,  unlike  the  auto¬ 
correlation  of  real  functions,  is  not  an  even  function  of  T.  Substituting 
from  (3)  for  G tx)  leads  to 


=  s(£). 

where  S('f)  is  the  Dirac  delta  function  and 

r  oo 

=  cU J  ^=r 


(5) 
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is  that  fraction  of  the  received  power  returned  by  scatterers  whose  velocity 
components  in  the  direction  of  the  transmitter  lie  in  t^U.  Since  a  frequency  f 
corresponds,  in  the  Doppler  sense,  to  a  velocity  component  toward 

the  transmitter,  the  result  (5)  can  be  expressed 


Ffr)*  SU). 


<M 


F(f)  is  defined  for  positive  and  negative  frequencies,  corresponding 
respectively  to  velocities  toward  and  away  from  the  transmitter.  Equation  (6) 
shows  that  the  power  spectrum  of  the  complex  amplitude,  henceforth  called 
the  Doppler  spectrum,  is  an  image  of  the  velocity  distribution  of  the 
scatterers,  weighted  according  to  the  cross  sections  of  the  scatterers  in  each 
velocity  interval.  This  simple  relationship  demonstrates  the  connection 
between  the  velocities  of  scatterers  and  the  radar  signal.  It  conveys  the 
same  information  as  equation  (3)  for  the  autocorrelation  function,  and  can 
be  used  as  the  basis  for  interpreting  signal  fluctuations. 

The  signal  sensed  at  the  antenna  is  not  the  complex  amplitude  eLir) 
but  is  the  quantity  £(f)  defined  by  (1).  It  is  of  interest,  therefore,  to  determine 
how  the  spectrum  of  Ertt)  ,  which  can  be  measured,  is  related  to  the 
Doppler  spectrum  of  the  complex  amplitude.  This  relation  can  be 

determined  by  first  noting  that 

_  V  tOut  ,  \  l  W0t 

aeC0«  e(t)e  +  e  (+)e  . 


Then  forming  the  product  EOOEtt+T)  ,  taking  time  averages,  and  trans¬ 
forming  to  the  frequency  domain  yield 
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Here  F 


■.«-r 


g(f)Eff-t-r)  3-27Tt^  T 


£(f) 


is  the  power  spectrum 


of  the  received  signal  E  (t)  •  The  following  normalizations  hold  for  the  spectra 
as  defined  above: 

f  Fm(fW  --  J «  I- 


Power  spectra  are  sometimes  defined  for  positive  values  of  frequency  only, 
because  the  detection  process  used  in  measuring  a  spectrum  makes  no  distinc¬ 
tion  between  positive  and  negative  frequencies  of  the  same  value.  The 
measured  spectrum  (f)  of  ECi)  can  therefore  be  defined  as 


Fjf)  (  z  *  F(-M)  t  F  f  so 


V  -  o;  f  <o. 


(7) 


with  the  normalization 


roc 

-  I- 

0 

The  spectrum  defined  by  (7)  is  available  for  measurement 

in  experiments.  The  Doppler  spectrum  F(f)  from  weather  targets  is 
restricted  to  an  interval  ,  in  the  region  of  f-Q  ,  which  is  very  small 

compared  to  the  carrier  F0  .  Consequently,  F(— 0  f°r  all  positive 
-f  ,  and  equation  (7)  shows  that  F^  (f )  is  simply  the  Doppler  spectrum 
displaced  along  the  frequency  axis  by  an  amount  ,  It  is  often  inconvenient 
to  spectrum-analyze  a  signal  containing  components  as  high  as  F0  .  In  this 
case  fc(f)  can  be  heterodyned  down  to  a  lower  effective  carrier  frequency, 
say  ■£'  ,  before  being  analyzed.  Then  the  measured  spectrum  is  given  by  (7) 
with  replaced  by  -f'  .  No  difficulty  arises  from  heterodyning  the  signal 
unless  -f'  is  chosen  smaller  than  £sf  ,  in  which  event  the  spectrum  will  have 
a  contribution  from  the  "folded"  component  .  Figure  1  illustrates 

this  folding  effect  on  the  measured  spectrum.  The  folding  is  most  severe 
when  the  received  signal  is  heterodyned  against  the  transmitted  signal. 
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In  this  case  0  and  any  asymmetries  in  the  Doppler  spectrum  are 
undetected.  If,  however,  the  velocity  distribution  and  hence  the  Doppler 
spectrum  were  even  functions,  Fv»  ff)  would  yet  be  sufficient  information 
for  describing  the  velocity  distribution. 

Transmitting  in  pulses  introduces  an  additional  consideration.  If  'ff 
is  the  pulse  repetition  frequency,  the  highest  possible  frequency  in  F&) 
is  ^  "fy  .  Consequently,  F(f)  will  be  an  image  of  the  velocity  distribution 
only  if  no  radial  velocity  components  greater  than  4"^^  are  present 
in  the  target. 

We  conclude  that  the  power  spectrum  of  the  heterodyned  signal  can  be 
used  to  calculate  the  radial  velocity  distribution  of  the  scatterers.  The  receiver 
must  be  linear  and  coherent  with  the  transmitter,  all  signal  processing  prior 
to  spectrum  analysis  must  be  linear,  and  the  final  effective  carrier  frequency 
must  be  chosen  large  enough  to  prevent  folding.  Then  the  relation  between 
the  measured  spectrum,  the  Doppler  spectrum,  and  the  velocity  distribution 
is 


f„w  =  fc f~r)  =  su-  ^)- 


(8) 


The  Doppler  spectrum  FvfJ  ,  as  we  have  defined  it,  is  an  intrinsic  property 
of  the  radar  signal  independent  of  the  way  the  signal  is  received  and  processed. 
When  using  Doppler  radar,  the  quantity  can  be  measured  and,  as 

shown  in  (8),  is  closely  related  to  the  Doppler  spectrum.  With  incoherent 
radar  it  is  impossible,  except  in  special  cases,  to  calculate  F(r)  (or 
correspondingly  the  velocity  distribution)  from  any  measured  quantities. 
Certain  properties  of  the  velocity  distribution  can  be  derived  from  incoherent 
measurements,  however,  and  it  is  convenient  to  describe  these  properties 
in  terms  of  the  Doppler  spectrum. 


A-8 


VC-1606-P-1 


4.  Signal  power 

The  power  1(f)  in  the  signal  described  by  (1)  is  proportional  to 

'  tf'it)  *  anc*  can  ke  measured  as  the  output  of  a  square-law  detector 
in  the  radar  receiver.  The  power  is  the  quantity  ordinarily  measured  by 
incoherent  weather  radars.  It  is  a  fluctuating  quantity  because  /Hf)  fluctuates; 
the  spectrum  of  its  fluctuations  is  connected  with  the  Doppler  spectrum  or 
5(i».)  ,  but  in  general  cannot  be  used  to  obtain  5^^)  because  the  phase  infor¬ 

mation  contained  in  Eii")  is  lost  in  the  process  of  forming  X  It)  .  We  will 
now  consider  how  the  loss  of  phase  information  restricts  the  data. 

Siegert  and  Fleisher  have  both  shown  that  the  autocorrelation  ^(t) 
of  the  power  fluctuations  is  related  to  the  autocorrelation  of  the  complex 
amplitudes  by 

^<V)  -  G-(-r)  «-*( t)  .  (9) 


Of  course  6(r)  is  generally  complex,  and  it  can  be  written  so  as  to  include 
the  real  function  <j.(f)  : 

£(t)=  fr)e 


where  f/tl)  is  also  real.  In  order  that  the  symmetric  property  of  £(t)  be 
satisfied,  namely  Gf'C'Z)  ,  it  is  necessary  and  sufficient  that  tyft)  be 

even  and  ifit)  be  odd. 

The  power  spectrum  SpCt)  of  power  fluctuations  is  obtained  by  trans¬ 
forming  the  autocorrelation  function.  According  to  (9),  equivalent  expressions 
for  the  spectrum  are 

Sp  Cf  U  Jf  «.*( r)  e~*m*"*d  x 
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Both  of  these  relationships  can  be  written  in  terms  of  power  spectra  rather 
than  correlation  functions  by  noting 

d+ 


and  by  defining 


These  substitutions  lead  to 

do  r  ^ 

sr  (fW  F(v^  F'c ^  ^  (10) 

as  equivalent  expressions  for  the  spectrum  of  power  fluctuations.  F(f)  is 
the  Doppler  spectrum  while  FT'fJ  is  an  even  function  of  frequency,  obtained 
by  transforming  .  The  function  F'Cf)  can  be  computed  from  measurements 

with  incoherent  radar.  It  is  of  interest  to  determine  under  what  circumstances 
F'Cf)  can  be  used  to  obtain  the  velocity  distribution. 

Suppose  the  velocity  distribution  is  an  even  function:  that  is,  5(u)  -  S(-u). 
Then  r(f)  =  F('f)  is  true  for  the  Doppler  spectrum.  This  implies  immedi¬ 
ately  that  G( t)  =-  .  We  already  know  that  ,  so  we  conclude 

Gi G-^Yt), 

or  G(t)  is  a  real  function.  Therefore,  from  (9), 

G(t)  =  j  (tJ. 

In  terms  of  the  power  spectra,  this  result  becomes 

F(F)=  F'(-f). 
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This  shows  that  the  Doppler  spectrum  is  the  same  as  the  Fourier  transform 
of  the  square  root  of  the  autocorrelation  of  power  fluctuations  when  the  velocity 
distribution  is  even.  The  velocity  distribution  can  therefore  be  determined  from 
incoherent  measurements  when  it  is  an  even  function. 

More  generally,  it  can  be  verified  that  symmetry  in  the  Doppler  spectrum 
(or  velocity  distribution),  rather  than  evenness,  establishes  a  ftrong  connection 
between  the  two  spectrum  functions.  If  F(f)  is  symmetrical  about  ^  ,  that 
is,  -  F(^— f)  ,  then 

This  result,  illustrated  in  Figure  2,  means  that  the  spectrum  obtainable  from 
incoherent  measurements  will  have  the  same  appearance  as  the  Doppler  spectrum 
when  the  velocity  distribution  is  symmetrical.  The  mean  velocity  however, 
which  corresponds  to  Doppler  frequency  ^  ,  is  lost  in  the  incoherent  measure¬ 
ment. 

5.  Level-crossings  of  the  detected  signal:  the  R-meter 

A  still  more  limited  amount  of  information  about  the  velocities  of  the 
scatterers  is  contained  in  the  rate  at  which  the  power,  or  some  function  thereof, 
crosses  an  arbitrary  level.  Rutkowski  and  Fleisher  (1955),  for  example, 
measured  the  rate  at  which  the  signal  from  a  linear  detector  crossed  its 
average  value,  and  interpreted  the  results  in  terms  of  velocities.  An  instru¬ 
ment  to  measure  this  rate,  called  the  R-meter,  has  been  used  by  several 
investigators  and  is  incorporated  on  some  Weather  Bureau  research  radars. 

Such  measurements  will  now  be  shown  to  be  related  to  certain  properties  of 
the  Doppler  spectrum. 

Rice  (1944)  proved  that  the  expected  rate  of  zero  crossings  for 
Gaussian  noise  is  uniquely  determined  by  the  spectrum  of  the  noise.  His 
analysis  can  be  modified  for  application  to  signals  that  are  not  Gaussian,  such 
as  detected  radar  signals.  The  power  returned  are  not  Gaussian,  such  as 
detected  radar  signals.  The  power  returned  from  random  scatterers  is 
distributed  as  the  sum  of  the  squares  of  two  equal  independent  Gaussian  random 
variables,  as  elaborated  upon  by  Lawson  and  Uhlenbeck  (1950).  Using  their 
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results  in  connection  with  the  technique  described  by  Rice  leads  to  the 
following  expression  for  the  expected  rate  L  at  which  the  power  crosses  its 
average  value: 


L- 


(Fn 
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[  /~r  M 
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where  is  the  power  spectrum  of  power  fluctuations.  Substituting  from 

(10)  gives 

\_-  -5—  f  F(»-rf)  F(>)d»]  , 


which  reduces  to 
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where  C£.  is  defined  by 
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In  terms  of  the  radial  component  It  of  velocity, 


Therefore,  the  mean  rate  at  which  the  received  power  crosses  its  average 
level  related  to  the  standard  deviation  of  the  velocity  distribution  by 


L- 


8jfTT 
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This  relationship  is  applicable  when  the  quantity  measured  is  proportional  to 
power,  i.e.  when  a  square  law  detector  is  used  in  the  receiver.  If  a  linear 
detector  is  used,  the  measured  quantity  is  then  proportional  to  [j ^ 
and  the  result  of  Rutkowski  and  Fleisher  (1955)  is  applicable:  the  mean  rate 
at  which  this  signal  crosses  its  average  is 

4-TT  -TT/4-  , 

—  e  ai  =  S-.74-  tfi, I ■ 
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Thus  the  proportionality  factor  which  connects  the  standard  deviation  of 
velocity  with  the  rate  of  average  level  crossings  varies  only  about  10  per  cen' 
as  the  power  law  of  the  detector  ranges  from  linear  to  square.  Small  deviations 
from  an  exact  linear  or  square  function  in  the  detector  are  apparently  of  little 
importance  in  determining  the  proportionality  factor. 
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ABSTRACT 


From  the  time-varying  Doppler  spectrum  of  certain  meteorological 
radar  targets  there  can  be  derived  the  energy  in  the  radial  component  of 
turbulence  at  target  altitude.  Two  properties  of  the  spectrum  must  be  re¬ 
corded  continuously:  the  mean  Doppler  velocity  and  the  variance  <T*  . 

Then  the  variance  %>  of  the  time-varying  mean  velocity  and  the  average 
value  (T of  the  variance  must  be  computed  from  the  record.  The  sum  of 
these  two  computed  quantities  is  proportional  to  the  total  energy  E  of  the 
radial  component  of  turbulence  per  unit  mass  of  air.  Specifically, 

*  £  [%)  +  ]  > 


where  Si)  is  the  three-dimensional  spectrum  of  the  radial  component  of 
turbulence. 
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1. 


Introduction 


The  radar  signal  backscattered  from  a  weather  target  or  a  cloud  of 
chaff  can  be  processed  *■'">  obtain  information  about  the  velocities  of  the  par¬ 
ticles  comprising  the  target.  Because  of  wind  shear,  turbulence,  or  unequal 
fall  velocities,  these  particles  are  characterized  by  a  velocity  spread  or 
distribution  rather  than  by  a  single  translational  velocity.  The  distribution 
of  the  relative  velocity  of  two  particles  chosen  at  random  can  be  obtained  by 
analyzing  the  fluctuations  in  power  received  by  conventional  noncoherent 
radars  (Stone  and  Fleisher,  1956).  For  viewing  conditions  in  which  the  effects 
of  wind  shear  and  unequal  fall  velocities  can  be  estimated,  the  spread  of  the 
distribution  of  relative  velocities  has  been  used  as  a  description  of  turbulence 
intensity,  usually  with  the  assumption  that  the  illuminated  particles  respond 
instantly  to  wind  changes  (Rutkowski  and  Fleisher,  1955). 

The  scales  of  turbulence  responsible  for  the  spread  of  the  velocity 
distribution  are,  in  a  rough  sense,  those  smaller  than  the  dimensions  of  the 
radar  target  volurr  - .  Larger  scales  of  turbulence  are  not  effectively  manifested 
by  spreading  in  the  relative -velocity  distribution,  but  rather  by  time  variations 
in  the  distribution  of  particle  velocities.  The  latter  distribution,  while  not 
obtainable  from  measurements  with  noncoherent  radars,  can  be  simply  de¬ 
rived  from  Doppler  (coherent)  radar  measurements. 

For  those  weather  targets  consisting  of  particles  that  respond  faithfully 
to  the  wind,  the  velocity  distribution  of  the  particles  can  be  interpreted  as 
that  of  the  wind.  Therefore,  since  it  measures  the  complete  distribution  of 
wind  speed  in  a  volume  in  space,  Doppler  radar  can  be  considered  as  a 
volume -sampling  anemometer.  The  purpose  of  this  paper  is  to  describe  a 
way  in  which  such  an  anemometer  can  be  used  to  measure  the  total  turbulent 
energy  of  the  wind  within  any  chosen  volume  of  the  atmosphere. 

2.  Relation  between  Doppler  Spectrum  and  Particle  Velocity  Distribution 

The  spectrum  of  the  fluctuating  signal  amplitude  received  by  a  Doppler 
radar  is  closely  related  to  the  distribution  of  radial  components  of  velocities 
of  the  scatterers  comprising  the  radar  target.  The  velocity  distribution  can 
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be  described  by  the  function  S(V),  where  S(V)dV  is  the  relative  number  of 
par'a*  les,  weighted  by  their  average  radar  cross  sections,  whose  radial 
velocity  components  lie  in  the  interval  dV.  This  function  is  defined  for  posi¬ 
tive  and  negative  V,  and  is  normalized  to  unit  area.  If  the  cross  sections 
and  velocities  of  the  scatterers  are  uncorrelated,  then  the  weighting  is  uni¬ 
form  and  S(V)  is  the  probability  density  of  radial  velocity. 

If  the  output  of  a  radar  receiver  which  is  linear,  and  coherent  with  the 
transmitter,  is  heterodyned  down  to  frequency  f'  and  then  spectrum-analyzed, 
a  measured  spectrum  Fm(f)  results.  This  spectrum  is  related  to  the  weighted 
velocity  distribution  S(V)  by 


FJf)~  +  s(~£  -4rj/  , 

=  0  otherwise, 

where  (f-f)  \  /  2  «  V  is  the  usual  relation  between  the  Doppler  shift  (f-f ' ), 

carrier  wavelength  \,  and  radial  velocity  V.  The  first  term  in  the  brackets 
would  be  expected  from  translation  of  the  spectrum  along  the  frequency  axis 
as  a  result  of  heterodyning.  The  second  term  arises  from  the  folding  effect 
which  occurs  when  the  translated  spectrum  does  not  vanish  at  f  =  0.  These 
relations  are  illustrated  in  the  Figure.  We  shall  always  choose  f'  large 
enough  to  prevent  folding,  in  which  case 

FrnC*)*  S(4r  -  , 

for  positive  f.  Thus  the  measured  Doppler  spectrum  is  simply  an  image 
of  the  velocity  distribution  when  frequencies  are  interpreted  in  the  Doppler 
sense. 

3.  Weather  Targets  and  the  Wind 

For  practical  expedience  in  measuring  the  wind  and  its  properties, 
we  would  like  S(V)  to  be  the  same  as  the  distribution  of  the  radial  component 
of  wind  in  the  radar- sampled  volume.  This  wind  has  a  mean  and  a  turbulent 
component,  both  of  which  would  ideally  be  represented  by  S(V).  Stackpole  (1961) 
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has  shown  that  raindrops  cannot  be  relied  upon  to  respond  well  to  gusts. 
Consequently,  the  measured  S(V)  from  a  target  of  raindrops  will  give  only 
biased  estimates  of  turbulence.  Cloud  droplets,  on  the  other  hand,  are 
ideal  wind  tracers,  but  may  not  always  backscatter  enough  energy  to  provide 
a  detectable  target  for  insensitive  radars.  Snowflakes  lie  somewhere  between 
raindrops  and  cloud  drops  in  terms  of  their  responsiveness  to  wind  and  they 
are  detectable  at  reasonable  ranges  by  radars  of  modest  capability.  It  has 
been  found  by  extending  Stackpole's  analysis  that  dry  snowflakes  respond 
well  enough  to  gusts  to  be  considered  adequate  tracers  for  both  the  mean 
and  the  turbulent  components  of  the  wind.  Additionally,  their  fall  velocities 
are  usually  small  -  less  than  1  m/sec  -  and  the  measured  S(V)  can  be 
adjusted  to  account  for  those  velocities.  Although  there  is  a  correlation 
between  the  cross  section  and  fall  velocity  of  a  snowflake,  the  range  of 
vertical  velocities  is  so  small  that  no  serious  error  is  introduced  by  neglec¬ 
ting  the  weighting  on  S(V)  due  to  cross  section  differences.  Thus  the  wind 
distribution  can  easily  be  derived  from  the  measured  velocity  distribution 
of  snowflakes. 

Let  us  denote  the  wind  distribution  in  the  radar-sampled  volume  by 
S'(V).  This  function  will  differ  noticeably  from  S(V)  only  for  vertical  viewing, 
when  the  full  component  of  fall  is  detected.  The  properties  of  S ' ( V)  which  are 
important  for  our  analysis  are  its  mean  value  and  the  spread  of  velocities 
about  that  value.  The  effect  of  small  scales  of  turbulence  is  to  spread  the 
distribution,  while  larger  scales  cause  the  location  of  the  distribution  on  the 
V-axis  to  change  with  time.  In  application  to  turbulence  measurements,  the 
analysis  is  simplified  if  there  are  no  small  scale  components  of  mean  wind 
shear  so  that  the  spread  of  S'(V)  is  entirely  due  to  turbulence;  otherwise  the 
contribution  to  spread  due  to  mean  wind  shear  must  be  subtracted  from  the 
total  spread.  The  subtraction  to  account  for  mean  shear  would  require  an 


Determination  of  the  fidelity  with  which  meteorological  targets  follow  the 
wind  requires  lengthy  computations  which  have  not  yet  been  carried  out  for 
typical  radar  chaff.  Presumably  they  could  be. 
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independent  and  precise  measurement  of  the  wind  in  the  sampled  volume,  and 
such  a  measurement  is  not  presently  possible.  Mean  wind  shear  does  not 
affect  measurements  for  vertical  viewing,  however,  and  the  small  spread 
introduced  into  S(V)  by  the  settling  of  snowflakes  can  easily  be  estimated  and 
subtracted  from  the  total  spread  to  obtain  S'(V).  Therefore  vertical  meas¬ 
urements  in  snow  constitute  at  least  one  situation  in  which  the  turbulent  wind 
distribution  can  be  obtained  from  the  radar  data. 

4.  The  Measurement  of  Turbulent  Energy 

As  indicated  above,  S'(V),  the  distribution  of  radial  components  of 
turbulent  velocity  within  the  sampled  volume,  is  readily  obtained  from  S(V)  for 
the  case  of  vertical  viewing.  The  time  variations  of  the  velocity  variance 
and  the  average  velocity  can  be  obtained  directly  from  S'(V)  using  analog 
techniques.  We  will  now  derive  a  fundamental  relation  between  these  quantities 
and  the  energy  of  the  observed  component  of  turbulence. 

The  (time)  average  value  of  the  variance  of  the  radial  component  of  flow 
velocity  V(r  )  in  the  sampled  volume  v  is  defined  by 


This  variance  is  identical  to  the  time  average  of  the  velocity  variance  of  S ' ( V) . 
Since  time  averages  commute  with  space  averages,  the  above  equation  can  be 
rewritten  as: 

*  £.(  v'(r)d«r  -  f  <t*i.  Vtf)V( K) 

•v  "vr  ^ 

*  V*(r)  ~ 

where  By  (r^  -  r^  )  is  the  three-dimensional  correlation  function  of  V(r  ), 


See  Tatarski,  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill,  1961, 
for  a  discussion  of  the  functions  used  in  describing  random  fields. 
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and  time  averages  have  been  replaced  by  ensemble  averages.  If  the  field  V(F) 
is  assumed  to  be  homogeneous,  the  three-dimensional  spectrum  of  V 

can  be  introduced,  giving 

(Ty  (*r)  =  fff$(2)[-  +('*’:*}]  <*%  >  (1) 

where  k  is  the  wave  vector  and  f(v;Tc)  it  a  function  determined  by  the  size  and 
shape  of  the  radar- sampled  volume.  For  any  volume,  this  function  equals  unity 
at  k  =  0  and  approaches  zero  as  k  increases.  Thus  <r^  («r)  is  more  sensitive 
to  large  wave  numbers  than  to  small  ones.  For  ordinary  radar  target-volumes, 
the  result  is,  roughly,  that  the  average  variance  of  the  velocity  distribution  is 
sensitive  to  scales  of  turbulence  smaller  than  a  characteristic  length  of  the 
volume,  and  relatively  insensitive  to  scales  larger  than  this  length. 

The  other  property  of  S'(V)  related  to  the  turbulent  energy  is  the  variance 
in  the  time  domain  of  the  average  value  ^V/  of  velocity  over  the  distribution. 

A  development  similar  to  that  above  shows  this  variance  to  be 


where  f(v;k  )  is  the  form  function  previously  described.  A  comparison  of  (1) 
and  (2)  shows  that  is  sensitive  to  the  form  function  f(v  ;"I? )  in  a  manner 

exactly  complementary  to  i.e.  »  the  quantity  defined  by  (2)  is  sensitive 

to  the  large  scales  of  turbulence,  and  relatively  insensitive  to  the  small  scales. 
As  already  indicated,  the  converse  is  true  for  <r^~(«r)  .  Indeed,  the  relative 
magnitudes  of  and  crv  («r)  indicate  the  partitioning  of  the  turbulent 

energy  between  turbulent  scales  larger  than  and  smaller  than  the  sampled 


volume. 


Adding  (1)  and  (2),  we  obtain  a  result  independent  of  the  sampled  volume. 


+  (T^(*r) 


The  term  on  the  right  is  equal  to  twice  the  energy  E  ,  per  unit  mass  of  air, 
in  the  measured  component  V  of  turbulence.  Consequently,  the  energy  is 
related  to  the  computed  quantities  by 
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The  total  turbulent  energy  E^,  (in  all  components)  can  be  obtained  only  with 
the  assumption  of  isotropy,  in  which  case 

=  i  ['%>("'>+  <t; VJ  7  • 
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Figure  I  RELATION  BETWEEN  VELOCITY  DISTRIBUTION  OF  SCATTERERS 
AND  MEASURED  SPECTRUM  AFTER  HETERODYNING 
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